Back

G3: Genes, Genomes, Genetics

◐ Oxford University Press (OUP)

Preprints posted in the last 30 days, ranked by how well they match G3: Genes, Genomes, Genetics's content profile, based on 252 papers previously published here. The average preprint has a 0.20% match score for this journal, so anything above that is already an above-average fit.

1
Large-scale pooled selection enables mapping of modifiers of RNAi gene knockdown effects in C. elegans

Zdraljevic, S.; Lee, A.; Walter-McNeill, L.; Marquez, H.; Webster, A. K.; Chitrakar, R.; Blount, J.; Baugh, L. R.; Bloom, J. S.; Kruglyak, L.

2026-10-11 genetics 10.64898/2026.10.08.757717 medRxiv
Top 0.1%
22.0%
Show abstract

Genetic interactions contribute to phenotypic variation, yet most studies characterize them in a small number of reference backgrounds rather than across natural populations. Studying genetic interactions in Caenorhabditis elegans is facilitated by RNAi by feeding, which enables gene-specific knockdowns across many genetic backgrounds without the need to construct individual mutant lines. However, individually screening hundreds of isolates across many RNAi conditions is laborious. Here, we developed a pooled platform that uses whole-genome sequencing to quantify the fitness of hundreds of genetically distinct C. elegans strains competing in a single flask. We leveraged this platform to evaluate RNAi responses in a pooled population of RNAi-sensitive C. elegans strains and identified two classes of QTL that influence RNAi responses. The first class influences responses to knockdowns of specific genes and is exemplified by a QTL that is specific to knockdown of mig-6. The second class influences general RNAi-response across multiple knockdowns. One QTL in this class localized to the right arm of chromosome III and was detected by both GWA and xQTL mapping, which suggests that a common large-effect variant underlies this QTL. We identified a T96K missense change in the dsRNA receptor SID-2 as the underlying causal variant. The 96K allele resides in the extracellular domain of SID-2 and confers increased RNAi sensitivity relative to the 96T allele. We hypothesize that the 96K allele enables greater uptake of ingested lumenal dsRNA through higher-affinity binding to the dsRNA phosphate backbone. The platform we developed enables systematic mapping of genetic interactions across natural populations at a previously impractical scale.

2
Fluorescent Ubiquitination-based Trbl Turnover Indicator (FUTTI): A tool to detect activity of the Tribbles pseudokinase in vivo

Nauman, C.; Shayestehpour, S.; O'Connor, S.; Fischer, Z.; Dobens, L.

2026-10-01 genetics 10.64898/2026.09.25.753906 medRxiv
Top 0.1%
19.2%
Show abstract

Targeted protein turnover of the transcription factor C/EBP is enhanced by the Tribbles (Trib) pseudokinase family of adaptor proteins, which bind to a conserved degron motif in C/EBP substrate and recruit a cognate COP1 E3 ligase to ubiquitinate C/EBP for degradation. We identified a conserved degron in the Drosophila C/EBP ortholog Slbo and fused this motif to a fluorescent protein (GFP or RFP) under the control of a constitutive promoter. Transgenic flies expressing this fluorescent protein (FP)::degron transgene, which we call the Fluorescent Ubiquitination-based Trbl Turnover Indicator (or FUTTI) sensor, showed reduced levels of fluorescence in response to increased levels of Drosophila Tribbles (Trbl) and conversely, increased levels when Trbl activity was reduced. Mutations in the Trbl conserved pseudokinase substrate binding pocket and in a divergent PEST-like motif in the C-terminal tail region disrupted Trbl activity in the FUTTI assay. These data indicate that FUTTI will be a valuable tool to identify sites of Trbl activity in vivo and the critical features of the protein required for function.

3
Domains of RNA-binding protein Ssd1 required for tolerating chromosome amplification in yeast

Gasch, A. P.; Jochem, A.; Ahrens, J. M.; Mayhew, E. M.

2026-09-28 molecular biology 10.64898/2026.09.25.754501 medRxiv
Top 0.1%
18.4%
Show abstract

Many RNA-binding proteins play pleiotropic roles in the cell, often through disparate domains that can be decoupled. Here we undertook a mutational analysis of yeast RNA-binding protein Ssd1 to better understand domains and processes required for tolerating extra chromosomes, a special class of aneuploidy. Ssd1 is important for cellular fitness when cells carry extra chromosomes, for reasons that are not clear. Ssd1 has multiple distinct domains, many of which are poorly characterized, including multiple RNA binding folds, nuclear import and export signals, an intrinsically disordered domain, an RNA polymerase II-interacting region, and many phosphorylation sites. We measured multiple phenotypes, including the ability to tolerate an extra copy of yeast chromosome XII as a representative aneuploid. Mutating RNA binding domains disrupts fitness in aneuploid, but not euploid, cells and ablates binding of a representative Ssd1 target. However, mutations affecting phase separation or interaction with the kinase Cbk1, which regulates Ssd1 granule formation during the cell cycle, did not affect aneuploidy tolerance. Evaluating correlations across these and other growth phenotypes informed on functional relationships among domains. We discuss implications for Ssd1 function and its role in aneuploidy tolerance in yeast.

4
Sorghum ANTHRACNOSE RESISTANCE GENE3 Is a Non-Coding RNA That Confers Fungal Resistance through Enhanced Cell Death

Mewa, D. B.; Adeyanju, A. O. A.; Guzha, A.; Okoye, P. C.; Liao, C.-J.; Seo, E.; Gribskov, M.; Lisch, D.; Mengiste, T.

2026-09-22 plant biology 10.64898/2026.09.18.752636 medRxiv
Top 0.1%
18.3%
Show abstract

The fungal pathogen Colletotrichum sublineola is a major constraint to sorghum production. The sorghum genotype IS18760 is resistant to multiple C. sublineola strains, but the genetic basis of this resistance is unknown. Bulk segregant analysis using whole-genome resequencing of a biparental mapping population identified a major resistance locus, designated ANTHRACNOSE RESISTANCE GENE LOCUS 3 (ARG3), associated with extensive cell death. Fine mapping delimited the ARG3 locus to a 30 kb interval containing three predicted protein-coding genes and one unannotated non-coding RNA gene. Among these, only the non-coding RNA gene showed consistent sequence polymorphisms between resistant and susceptible lines and elevated pathogen-induced expression in the resistant parent, suggesting that it underlies ARG3-mediated resistance. Regions syntenic to ARG3 across grass species contain conserved non-coding sequences enriched for H3K27me3 and H2A.Z chromatin marks, suggesting a conserved regulatory function. The resistance allele is absent from a 111-line sorghum pangenome. The resistant sorghum line harboring ARG3 and the corresponding avirulent pathogen strain originate from the same geographic region, suggesting co-evolution between host resistance and a prevalent pathogen strain. Phylogenetic analysis suggests that loss of resistance in susceptible lines may relate to insertion of a 12.4 kb repetitive sequence. Importantly, silencing of the ARG3 transcript in the resistant parental line abrogated resistance to C. sublineola, confirming that ARG3 non-coding RNA gene underlies the fungal resistance mediated by the ARG3 locus in IS18760. These lines of evidence show that ARG3 is an important target for resistance breeding and provide new insight into non-coding RNA-mediated disease resistance.

5
ABC transporters support Drosophila intestinal stem and progenitor cells during cytotoxic challenge

Dayton, H.; Russo, J.; Kolbert, K. E.; Williamson, O.; Balciunaite, A.; DSouza, E.; Becker, K.; Hosage, E.; Issa, M.; Lu, V.; Viswanatha, R.; Kondo, S.; Perrimon, N.; Markstein, M. M.

2026-10-12 genetics 10.64898/2026.10.05.756803 medRxiv
Top 0.3%
11.9%
Show abstract

Adult stem cells can survive diverse insults, from ionizing radiation to toxic chemicals. In vertebrates, ABC transporter expression and dye-efflux capacity have been associated with stem-cell protection and multidrug resistance. How compound accumulation differs between adult stem cells and their differentiated daughters within an intact tissue remains poorly understood. Here, we show that adult Drosophila intestinal stem/progenitor cells maintain a distinct low-accumulation state for membrane-permeable compounds and that ABC transporters contribute to this phenotype and to the response to cytotoxic stress. Across nine membrane-permeable probes, stem/progenitor cells accumulated eight at lower levels than differentiated enterocytes. We selected tetramethylrhodamine methyl ester (TMRM), a membrane-permeable fluorescent probe that showed a strong difference in accumulation between stem/progenitor cells and enterocytes, as a quantitative readout for an RNAi screen of all 56 Drosophila ABC-family genes. The screen identified seven candidate transporters whose depletion increased the stem/progenitor-to-enterocyte TMRM fluorescence ratio. Further characterization of one of them, ABCB7, showed that ABCB7-depleted guts contained fewer mitotic cells during cytotoxic drug exposure. These results implicate ABC transporters in lower compound accumulation in Drosophila intestinal stem/progenitor cells and support a role for ABC transporters in their response to cytotoxic challenge.

6
Inbred and outcrossed Drosophila embryos diverge in stress defense and developmental maintenance during heat stress

Natua, S.; Biswas, D.; Rashid, F.; Kalsotra, A.; Sokac, A. M.

2026-10-06 genetics 10.64898/2026.10.01.754988 medRxiv
Top 0.3%
11.8%
Show abstract

Early embryos must mount stress defenses without compromising the developmental programs that drive morphogenesis, but whether genetic background alters this balance remains unclear. We compared cellularizing Drosophila melanogaster embryos from inbred and outcrossed backgrounds under control conditions and after acute or chronic heat-stress. Both backgrounds developed successfully and showed modest transcriptomic differences at 18{degrees}C. Heat stress at 32{degrees}C, however, sharply separated their phenotypic and transcriptional responses. Inbred embryos showed poor larval hatching after acute and chronic heat, and this deficit persisted through adult eclosion. Their transcriptomes exhibited a greater stress-response burden, with enrichment of proteostasis, DNA repair, mitochondrial maintenance, RNA processing, and metabolic remodeling programs. Outcrossed embryos retained substantially greater developmental survival and higher abundance of transcripts associated with development, adhesion, epithelial organization, signaling, transport, sensory responses, and morphogenesis. Background-dependent divergence was greatest during acute heat. Under chronic heat, background-dependent differences were fewer, with most overlapping those seen during acute heat and a smaller set detected only under chronic heat. Transcript-origin analysis further showed that acute outcrossed-biased genes were predominantly maternally derived and showed no enrichment for purely zygotic transcripts, whereas chronic outcrossed-biased genes contained a larger zygotic component. In summary, the less viable inbred background mounted the broader molecular stress signature, whereas the more viable outcrossed background retained greater abundance of development-associated transcripts. These findings identify genetic background as a major determinant of how embryonic heat responses are distributed between molecular stress defense and developmental gene expression.

7
Selection on cuticle color produces a genome-wide response that does not localize to known major genes in Drosophila

Ruckman, S. N.; Long, A. D.

2026-10-05 genetics 10.64898/2026.09.29.755441 medRxiv
Top 0.3%
11.8%
Show abstract

Complex traits are often deemed highly polygenic when a mapping study fails to localize the loci underlying them, but a result that fails to identify major genes can also reflect a lack of power. Cuticle color in Drosophila offers a natural test of this distinction, since genes of large effect segregate in natural populations. We selected for dark and light thoracic trident color for 16 generations in two populations of D. melanogaster and two of D. simulans, then used pooled sequencing to scan each genome for allele-frequency differentiation between selected and control lines. Selection produced a strong, genome-wide shift in allele frequencies in every population, exceeding the expectation under drift. Simulations demonstrate that our design has low power to detect major-effect genes, which may explain why we do not observe peaks localizing to known color genes. We further show that a highly significant peak in each genome, on chromosome 3L, is an artifact of chorion-gene amplification in the pooled adult females. A source of false signal in pool-seq studies that has not been reported before. Cosmopolitan inversions shifted in frequency under selection in the one population that carried them, but direct testing showed they did not inflate genome-wide noise, contradicting a common assumption in the field. Together, these results demonstrate that a diffuse, polygenic-type response can arise from a trait with known large-effect genes, and provide a quantitative account of why an evolve-and-resequence design can fail to localize even a well-characterized genetic architecture.

8
HapSolo2: GPU-accelerated optimization for removing secondary haplotigs from diploid genome assemblies

Yin, J.; Agrawal, M.; Goel, P.; Solares, E. A.

2026-10-08 bioinformatics 10.64898/2026.10.02.756147 medRxiv
Top 0.3%
10.8%
Show abstract

BackgroundSecondary haplotigs retained in a primary assembly inflate genome size and duplicate gene content, and their removal remains a necessary curation step for heterozygous diploid genomes. HapSolo (Solares et al. 2021) treats this step as an optimization problem: it searches three pairwise alignment thresholds (percent identity, query coverage, and the ratio of aligned query to aligned reference length) to minimize a cost defined over single-copy, duplicated, and missing conserved orthologs. The original implementation depended on a separate BUSCO v3 run for every contig, on Python 2.7 and pandas in its inner loop, on a single randomized forward-walk optimizer, and on OrthoDB v9 lineage datasets that are no longer distributed. ResultsHere we report HapSolo2, a reimplementation that replaces the per-contig BUSCO step with miniprot-based ortholog classification against OrthoDB v9, v10, or v12.2 lineage datasets, adds steepest descent and simulated annealing to the original random forward walk, and evaluates all walkers in batched CuPy array operations on a single GPU. On an avocado assembly of 5,122 contigs scored against 1,440 orthologs, 100,000 iterations of random walk and steepest descent ran 5.6- and 7.6-fold faster on one NVIDIA RTX A6000 than on 10 CPU processes. Across 45 runs on four diploid assemblies (Anopheles funestus, Persea americana, Vitis vinifera, and Amblyraja radi-ata; 212 Mb to 3.2 Gb) and three OrthoDB releases, the three optimizers reached best costs within 1.8% of one another. Independent assessment with BUSCO v6.1 showed that HapSolo2 reduced duplicated orthologs from 52.0% to 5.6% in P. americana, from 17.9% to 10.4% in V. vinifera, and from 7.8% to 2.8% in A. funestus, with completeness decreasing by at most 1.3 percentage points. On the avocado assembly, HapSolo2 reduced the cost reached by the original HapSolo from 0.2200 to 0.1036 and raised contig N50 from 3,366 kb to 4,608 kb. ConclusionsHapSolo2 makes ortholog-guided haplotig removal practical in minutes on a single workstation, independent of a BUSCO installation, and reproducible across OrthoDB releases. It is distributed as an installable Python package with a unified command-line pipeline, an assembly statistics and plotting tool, an Apptainer container, and archived copies of the legacy ortholog datasets.

9
Genome assembly of Eucalyptus recurva provides insights into inbreeding and conservation priorities in Australia's rarest Eucalyptus

McMaster, E. S.; Jones, A.

2026-09-25 genomics 10.64898/2026.09.21.753279 medRxiv
Top 0.3%
10.6%
Show abstract

Eucalyptus recurva (Mongarlowe Mallee) is Critically Endangered, with only six known adult individuals persisting across two sites in the Southern Tablelands of New South Wales, Australia. Its extreme rarity, uniquely long lifespan, and limited reproductive output make it a priority for conservation genomics. Here we report the first genome assembly of E. recurva: a haplotype-resolved, gapless, telomere-to-telomere (T2T) assembly produced from Oxford Nanopore Technologies (ONT) long-read sequencing. Both haplotypes span 11 chromosomes (consistent with the conserved Eucalyptus karyotype of 2n = 22), with assembly sizes of 521.1 Mb (Hap 1) and 501.5 Mb (Hap 2), BUSCO completeness >99.6%, and quality values of >QV 62. Comparative analyses place E. recurva within section Maidenaria as sister to Eucalyptus viminalis and reveal high synteny between the two species. Inter-haplotype comparison identified 3.6 million SNPs and modest structural variation, suggesting that, despite extreme demographic bottlenecking, E. recurva retains meaningful genomic heterozygosity. We also characterise the chloroplast genome, which exhibits heteroplasmy, and report an apparently bipartite mitochondrial genome. This reference genome provides an essential resource for conservation management, population genetics, and the study of eucalypt genome evolution.

10
Pooled progeny sequencing effectively estimates polyandry in species with large families

Farmer, T.; Bentz, A.; Macdonald, S. J.; Unckless, R.; Kelly, J.

2026-10-10 genomics 10.64898/2026.10.02.756305 medRxiv
Top 0.4%
9.8%
Show abstract

Polyandry, multiple mating by females, both benefits females by increases genetic variation within families and provides an opportunity for sexual selection through sperm competition and cryptic female choice. Estimating polyandry in species with large clutches is challenging because genotyping all progeny individually can be impractical. We present a likelihood-based framework for estimating the number of sires and their relative contributions from the maternal genotype and pooled sequencing of entire progeny pools. The method evaluates the probability of observed pooled offspring read counts under alternative paternal genotype mixtures, accounting for Mendelian segregation and sampling variation from finite sequencing depth. We implement this framework for both biallelic SNPs and short-read mini-haplotype markers constructed from closely linked SNPs, which provide multi-allelic information comparable to traditional microsatellite markers. Simulations show that this approach accurately recovers sire number and sire contributions across a range of conditions. We then applied the method to 275 wild-caught Drosophila melanogaster females and their offspring from a Kansas population using 3,053 biallelic SNPs and 239 four-allele mini-haplotype markers. Estimates from the two marker sets were highly concordant, and combined analyses indicate that 57% of families were sired by a single male. These results demonstrate that pooled progeny sequencing provides an efficient genomic approach for estimating polyandry in species with large families.

11
A genetic screen to identify novel interactions of the phosphatidylinositol phosphate signaling pathway in Drosophila spermatogenesis

Rossini, S.; Nasirzadeh, N.; Icmat, M. S.; Chu, F. H.-K.; Rana, M.; Micovic, K.; Singh, K.; Zhuo, J.; Shao, L.; Ohri, A.; Huang, K.; Lightwala, M.; Nagree, M. S.; Jiwani, T.; Bellamkonda, K.; Prasad, P.; Antonyan, L.; Lands, H.; Sesso, G.; Liu, L.; Liu, N. F.; Philson, A.; Sheykhkarimli, D.; Krithivasan, P.; Papadopoli, A.; Fabian, L.; Brill, J. A.

2026-10-01 genetics 10.64898/2026.09.25.754494 medRxiv
Top 0.4%
9.7%
Show abstract

Phosphatidylinositol (PI) phosphates (PIPs) are membrane lipids with crucial roles in intracellular signaling, membrane trafficking and cytoskeletal organization. To perturb PIP function in Drosophila melanogaster sperm development, we previously expressed the Salmonella PIP phosphatase SigD (also known as SopB) under control of a spermatocyte- specific promoter. SigD expression results in reduced levels of plasma membrane PI 4,5-bisphosphate [PI(4,5)P2; henceforth PIP2] and male sterility. Transgenic lines expressing high levels of SigD (SigD-high) exhibit defects in spermatocyte cytokinesis, ciliary transition zone maturation and spermatid elongation, whereas lines expressing low levels of SigD (SigD-low) have defects in nuclear shaping and spermatid cell polarization. To further investigate the role of PIPs during Drosophila sperm development, we conducted a sensitized genetic screen to identify suppressors and enhancers of SigD-low phenotypes. From a screen of 760 chromosomal deficiencies spanning 83.6% of chromosomes 2 and 3, we identified 80 interacting regions that were further narrowed down to 40 individual suppressors. Among the strongest interactors were genes involved in microtubule dynamics, nucleocytoplasmic transport and membrane trafficking. The novel genetic interactions identified by our screen pave the way for further investigation into molecular mechanisms underlying the affected aspects of sperm development.

12
Gap-free telomere-to-telomere haplotypes assembly of the hybrid grouper 'Yushuban' (Cromileptes altivelis female x Epinephelus fuscoguttatus male)

Cai, Y.; Wang, C.; Li, M.; Yang, Y.; Cai, C.; Lu, S.; Chen, S.

2026-09-25 genomics 10.64898/2026.09.20.752973 medRxiv
Top 0.4%
9.7%
Show abstract

The hybrid grouper 'Yushuban' (Cromileptes altivelis [female] x Epinephelus fuscoguttatus [male]) has shown promising performance in aquaculture trials, yet no complete reference genome has been available for this hybrid. We generated two gapless, telomere-to-telomere (T2T) haplotype-resolved genome assemblies of a single individual by integrating short-read, PacBio HiFi, Oxford Nanopore ultralong-read and Hi-C data. The maternal (Yushu_Calti_1.0) and paternal (Yushu_Efusc_1.0) haplotypes are 1,077.33 Mb and 1,083.36 Mb in length, with N50 values of 45.46 Mb and 45.55 Mb, respectively. Across all platforms, over 99.4% of reads mapped to each assembly, and Merqury-based quality values reached 57.7413 and 59.2068 (PacBio). BUSCO recovered 99.78% and 99.75% of complete core genes. We predicted 27,248 and 27,566 protein-coding genes, of which 98.52% and 99.23% were functionally annotated, and annotated telomeric motifs and repetitive elements (46.12% and 47.18% of each genome). These assemblies fill a gap in grouper genomic resources and provide a foundation for genetic dissection, parental genome comparison and molecular breeding of this hybrid.

13
Optimizing genome assembly, chromosome synteny, and genetic variant discovery from Oxford Nanopore sequences of Dermacentor reticulatus ticks

Dillon, K. C.; Sprong, H.; Ronai, I.; Choudhary, S.; Grant, C.; de Paula Baptista, R.; Ray, D. A.; Glenn, T. C.

2026-09-24 genomics 10.64898/2026.09.18.752137 medRxiv
Top 0.4%
9.6%
Show abstract

Generating high-quality genome assemblies for small animals with large genomes is complex due to their small body size, DNA contamination, and repetitive elements. Ticks exemplify these complexities, while also being a global health threat to humans, domestic animals, and wildlife. Advances in long-read sequencing platforms now make it feasible to obtain large amounts of raw sequence data from individual specimens, but challenges remain. Key genome workflow challenges include error correction, assembly, transposable element annotation, and chromosome assignment and synteny. Here we examine three individual Dermacentor reticulatus ticks using deep Oxford Nanopore sequencing. Comparing and contrasting bioinformatic tools for raw read and assembly manipulation allows us to identify the parameters that provide a high-quality haploid genome assembly among 38 assemblies. We find Dorado error corrected raw read data from approximately two flow cells are needed, but limiting assembly coverage to 40x produces the highest quality assemblies. Then we examined the effects of optimal workflow parameters on downstream analyses of gene synteny and manual transposable element annotation of the three tick assemblies in comparison to an independent assembly of a fourth D. reticulatus individual from a different country. Gene synteny analysis allows chromosome assignment of scaffolds and transposable element identification was improved markedly with manual curation. Finally, we compared the genetic variation of D. reticulatus across two populations and found similar genetic diversity. Our study provides a clear workflow to obtain high-quality assemblies from Oxford Nanopore sequences and genetic characterization of a species with a large, complex, and repetitive genome.

14
Exploring roles for essential proteins in yeast filamentous growth identifies the WASP homolog Las17 as a regulator of the Cdc42-dependent fMAPK pathway

Pujari, A. N.; Priyadarshini, A.; Li, Z.; Darekar, A. S.; Williams, D.; Climie, D.; Bahr, S.; Friesen, H.; Li, M.; Oken, J. M.; Romito, R. M.; Andrews, B. J.; Boone, C.; Cullen, P. J.

2026-09-24 genetics 10.64898/2026.09.21.751561 medRxiv
Top 0.4%
9.3%
Show abstract

Cell differentiation into distinct cell types generates functional specialization in eukaryotic organisms. Many fungal species undergo filamentous growth, where cells differentiate into elongated and adhesive filaments capable of expansion and invasion into new environments. The budding yeast Saccharomyces cerevisiae also undergoes filamentous growth, and the regulatory pathways that control the response have been well characterized; however, the roles essential proteins play in this response have not been systematically explored. To address this gap in understanding, we constructed a collection of 332 temperature-sensitive (ts) alleles in 320 essential genes in a strain background that undergoes filamentous growth (called {Sigma}1278b or Sigma). The ts Sigma collection showed differences in temperature sensitivity compared to a ts laboratory strain collection, revealing unexpected phenotypic diversity in essential alleles across populations of individuals. Screening the ts Sigma collection for phenotypes related to filamentous growth uncovered new phenotypes for >35% of essential alleles. New regulators of the Mitogen-Activated Protein Kinase (MAPK) pathway that regulates filamentous growth (fMAPK) were identified, including Las17, a homolog of Wiskott-Aldrich Syndrome Protein (WASP) in humans. Las17 regulated the fMAPK pathway by promoting delivery of the sensor protein, Sho1, and Rho GTPase Cdc42 to the plasma membrane. Las17 also functioned as a hub coordinating separate parallel aspects of the filamentation response. The widespread roles for essential proteins in regulating a eukaryotic differentiation response suggest broader roles essential proteins than is currently appreciated.

15
Ethanol stress induces siRNA-mediated epimutations and phenotypic variability in fission yeast

Villalba de la Pena, M.; Sarkies, P.

2026-10-09 genetics 10.64898/2026.10.04.756277 medRxiv
Top 0.5%
8.1%
Show abstract

Epigenetic information is required for the maintenance of gene expression states through cell division. Some epigenetic information can also persist through generations, raising the possibility that epigenetic differences may arise in populations and contribute to evolutionary processes such as drift and natural selection. Evaluating this requires characterising the properties of such epimutations, so that their rate, stability and genome-wide distribution can be assessed. Recent experiments in many organisms have demonstrated that epimutations occur, but have limited stability relative to DNA sequence mutations. However, the extent to which epimutations contribute to evolution is still unclear, in part because it is not clear whether epimutations respond dynamically to environmental conditions. Here, we use the fission yeast S. pombe, in which evolutionarily relevant epimutations have been observed, but not yet systematically characterised, to test whether a physiologically relevant stress, ethanol exposure, leads to changes in the properties of epimutations. Using laboratory evolution experiments with repeated population bottlenecks to minimise selection, we show that epimutations associated with changes in small interfering (si)RNAs occur at a greater rate in ethanol. Ethanol-induced epimutations are enriched at transposable element loci and correlate with changes in transposon expression. Ethanol exposure also increases gene expression noise genome-wide and variability in growth rate; importantly, both these increases are suppressed in mutants lacking siRNAs. Together, these results demonstrate that stress can alter the generation of heritable epigenetic information, and suggest that siRNAs may have a role in promoting phenotypic variability in response to environmental stress.

16
Identifying, phasing, and structurally annotating sex chromosomes for genome assemblies using CBS-tools

Whitt, L.; Akozbek, L.; Bentz, P. C.; Armstrong, E.; Harkess, A.; Carey, S. B.

2026-09-29 genomics 10.64898/2026.09.23.753885 medRxiv
Top 0.5%
7.9%
Show abstract

A complete reference genome for species with chromosomally-determined separate sexes should contain scaffolds for all sex chromosome homologs. However, sex chromosomes present distinct computational challenges compared to autosomes. Here we present a k-mer based analysis that utilizes whole-genome sequencing of a few sex-identified isolates: Cytogenetics-By-Sequencing (CBS) tools. Unlike other approaches that typically address one aspect of the sex chromosomes, CBS-tools strives to guide users from the discovery of the heterogametic sex through identifying the sex-determination region (SDR). The core of CBS-tools is automated quantification of sex-specific k-mers in order to predict the heterogametic sex. Using publicly-available datasets, CBS-tools correctly identified the known sex-system of the 31 species tested. Additionally, we used these k-mers to verify and correct phasing of sex chromosomes between haplotypes in species representing different sex-systems. Finally, we used these k-mers to delimit the SDR boundary using an interactive web platform. CBS-tools was developed with previously unexplored sex chromosome systems in mind, but is also suitable for well-examined sex chromosome pairs.

17
The Genetic Landscape of Structural Variants in Aspergillus flavus using ONT and PacBio Sequencing

Woods, A.; Lohmar, J. M.; Malik, A.; Lustig, A. J.

2026-10-11 genomics 10.64898/2026.10.05.756768 medRxiv
Top 0.6%
7.6%
Show abstract

Genetic rearrangement in the absence of selection is poorly understood. Long-read, high-depth sequencing has enabled analysis of the genetic landscape of spontaneous structural variants. We have applied PacBio Revio sequencing to identify the major classes of genetic alterations during long-term vegetative growth in Aspergillus flavus. These classes include a) precise deletions frequently associated with the flanking region of homology, b) clusters of short deletions that are heavily enriched within centromeric and subtelomeric regions that are typically 94% A+T, and c) deletions and expansions within a subset of homopolymer tracts. We find 70 homopolymer tracts longer than 44 bp, dispersed throughout the genome in noncoding regions and unassociated with transposable elements. A broader analysis of third-generation sequenced datasets from fungi found within the Phylum Ascomycota showed a high but variable prevalence of homopolymeric tracts across species. The adventitious identification of a putative flbA mutation, responsible for an altered conidial phenotype, by the Deep Variant AI system also demonstrates the power of this approach. These studies show that high-depth, long-read sequencing can probe SNPs, indels, homopolymeric tracts, and large-scale rearrangements in an unbiased fashion, and demonstrate that short deletions and homopolymer tract expansions arise spontaneously in A. flavus and other filamentous fungi.

18
A chromosome-level genome assembly of the chilli thrips Scirtothrips dorsalis Hood (Thysanoptera: Thripidae)

Busuulwa, A.; Liesenfelt, T.; Mongue, A. J.; Lahiri, S.

2026-10-10 genomics 10.64898/2026.10.08.757792 medRxiv
Top 0.6%
7.2%
Show abstract

Thrips are known for being pests of many agricultural crops, as well as vectors of plant viruses. Yet despite their economic importance, genomic resources for this group remain limited, and no chromosome-level genome has been available for the genus Scirtothrips. One such species is the chilli thrips, Scirtothrips dorsalis Hood (Thysanoptera: Thripidae). This invasive and highly polyphagous pest is established across much of the world and transmits several tospoviruses. Genomic resources would provide a new means to understand the traits that make thrips difficult to manage. Here, we present a chromosome-level genome assembly for S. dorsalis. The assembly is 266.5 Mb, with 96.9% of the sequence anchored in 16 chromosome-scale scaffolds. It contains 20,557 genes, and repetitive elements account for 28.8% of the sequence, which is comparable to those of other sequenced Thripidae. A species tree based on 3,019 single copy orthologs placed S. dorsalis as sister to Scirtothrips hansoni, with Scirtothrips most closely related to Frankliniella, Megalurothrips, and Odontothrips. The mitochondrial genome is a single circular molecule of 15,207 bp in which the sequence of the published mini circle is integrated, in contrast to the bipartite mitogenome previously reported for the same cryptic species.

19
Haplotype-phased, chromosome-level, annotated genome assembly of the soybean rust pathogen Phakopsora pachyrhizi

Greatens, N. F.; Elmore, M. G.; Ramachandran, S. R.; Wright, A.; Tallon, L. J.; Whitham, S. A.; Koch Bach, R. A.; Pedley, K. F.

2026-10-05 genomics 10.64898/2026.09.29.755439 medRxiv
Top 0.6%
7.1%
Show abstract

Phakopsora pachyrhizi, the causal agent of soybean rust, is among the most aggressive and economically impactful pathogens affecting soybean (Glycine max) production worldwide. As a dikaryotic rust fungus, P. pachyrhizi maintains its genetic information as two discrete haplotypes that are sequestered within separate nuclei for most of its life cycle as it reproduces asexually via the production of infectious urediniospores. Previous genomes of three P. pachyrhizi isolates collected in South America revealed a large genome with a high transposable element content, measuring up to 1.25 Gb and 93%, respectively. Although these efforts provided valuable foundational resources, the complexity of the genome resulted in highly fragmented assemblies. Here, we generated a highly contiguous, substantially improved, annotated genome resource from the P. pachyrhizi isolate LA04-1, collected during the initial outbreak of soybean rust in the continental United States in 2004. The assembly is haplotype-phased, with each haplotype assembled into 18 chromosomes totaling 572 Mb and 569 Mb in size. The isolate contains unique mating-type genes and does not share haplotype genomes with any available isolates with whole-genome sequencing datasets. This resource will be valuable to future investigations into the underlying mechanisms of virulence in this important plant pathogen.

20
Neuropeptide T2A-GAL4 knock-ins illustrate practical considerations for interpreting GAL4 reporter patterns in Drosophila

Jinnai, K.; Ozawa, K.; Yajima, K.; Tanimoto, H.; Kondo, S.

2026-09-22 genetics 10.64898/2026.09.16.750837 medRxiv
Top 0.6%
7.0%
Show abstract

Transgenic GAL4 lines are widely used in Drosophila to label and manipulate specific cell types, but the resulting reporter expression patterns do not always represent ongoing driver activity. We generated a collection of neuropeptide T2A-GAL4 knock-ins and characterized their expression by native-fluorescence imaging using UAS-CD8::GFP and a Brp-miRFP680 knock-in, eliminating the need for immunostaining. During characterization of these lines, we observed cases in which adult GFP patterns did not match the on-going GAL4 activity, as well as cases in which GAL4 produced adverse cellular effects. Using temperature-sensitive GAL80, we found that adult GFP expression in certain neurons driven by CCHa2-T2A-GAL4 and Dsk-T2A-GAL4 depended on GAL4 activity before eclosion and was not reproduced by adult-only induction. Several GAL4 lines also showed dominant developmental phenotypes. Heterozygous Burs-T2A-GAL4 animals, for example, displayed wing-expansion defects that were fully suppressed by GAL80. Severe wing defects were not accompanied by a reduction in the number of reporter-labeled neurons before eclosion, but the morphology of these neurons were affected. These observations illustrate how developmental driver history and GAL4-associated cellular effects can influence reporter patterns in the adult nervous system.